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GYROTRONS  AND  FREE  ELECTRON  LASERS 
.  FOR  ATMOSPHERIC  SENSING 


I.  Introduction 

Gyrotrons  have  been  developed  over  the  last  decade  and  a  half 
or  so  for,  and  to  a  large  extent,  by  plasma  physicists.  Free  electron 
lasers  (FEL's)  were  also  developed  in  part  by  plasma  physicists. 
Gyrotrons  have  been  developed  mostly  for  heating  of  fusion  plasmas, 
and  free  electron  lasers,  as  a  powerful  source  of  infrared  (IR) 
radiation  where  none  have  previously  existed.  Both  of  them  are 
generally  thought  of  being  situated  in  large,  fixed  facilities.  In  fact, 
FEL’s  often  are  the  fixed  facility,  as  they  require  a  high  energy 
electron  accelerator  and  a  heavily  shielded  enclosure.  However 
gyrotrons  are  compact  enough  that  they  could  be  fielded  in  for 
instance  a  small  truck.  While  IRFEL’s  cannot  currently  be  fielded, 
there  are  several  possible  development  strategies  which  could  lead 
to  a  Beldable  source  in  the  not  too  distant  future. 

The  field  of  atmospheric  sensing  has  benefited  enormously 
over  recent  decades  from  the  use  of  radar  and  laser  probes.  Radars 
are  generally  at  low  frequency,  typically  X-Band  or  below,  and 
tunable  lasers  are  usually  optical  or  shorter  wavelength.  (Fixed 
frequency  or  line  tunable  lasers  are  available  at  a  few  wavelengths 
in  the  infrared.)  The  development  of  fixed  frequency  gyrotrons 
opens  up  the  possibility  of  cloud  radars.  Millimeter  wave  radiation 
at  94  GHz  has  the  unique  capability  of  scattering  weakly  from  cloud 
aerosol  droplets,  which  really  define  the  cloud,  while  at  the  same 
time  propagating  easily  through  the  cloud.  It  could  for  instance 
sense  multiple  cloud  layers,  whereas  a  lidar  would  be  stopped  by  the 
first  cloud  not  visibly  transparent.  Tunable  gyrotrons  open  up  new 
possibilities  of  absorption  measurements  over  long  path  lengths. 
Detection  of  trace  elements  at  sea  level  and  in  the  upper  atmosphere 
is  also  possible.  Furthermore,  high  power  gyrotrons,  in  bistatic  radar 
configurations,  open  up  additional  possibilities  for  remote  sensing  of 
clear  air  turbulence.  Tunable,  powerful,  infrared  free  electron  lasers 
could  allow  for  the  detection  of  new  classes  of  atmospheric 
impurities  by  using  them  as  differential  absorption  lidars  at  much 
greater  range. 


Manuicr^  qiproved  Januuy  13,  1994. 


II.  Gyrotrons,  Free  Electron  Lasers  and  Atmospheric  Propagation 

High  power  gyrotrons  and  Free  Electron  Lasers  (FEL's)  have 
been  developed  in  Ae  United  States  and  around  the  world  for  the 
last  15  years.  In  each  case,  the  average  radiated  power  levels  now 
available,  respectively  in  the  millimeter  wave  and  infrared  regimes, 
are  far  in  excess  of  what  had  been  previous  capability.  At  this  point, 
gyrotrons  are  the  most  powerful  sources  of  radiation  at  millimeter 
wavelength  and  FEL's  have  achieved  high  power  in  the  infrared. 
While  this  radiation  is  generally  absorbed  in  the  atmosphere,  there 
are  propagation  windows.  In  Fig.(l)  is  shown  the  absorbtivity  of  the 
atmosphere  in  the  millimeter  wavelength^,  and  in  Fig.  (2a)  is  shown 
the  absorbtivity  in  the  infrared^.  The  white  areas  correspond  to 
propagation.  (These  curves  are  just  rough  rules  of  thumb,  the  actual 
propagation  depends  strongly  on  conditions.  Humidity  greatly 
reduces  propagation  in  the  millimeter  wave  regime,  and  infrared 
radiation  cannot  propagate  through  clouds  or  fog.)  Also  shown  in 

Fig.(2b)  are  the  approximate  capabilities  of  various  tunable  lasers. 
(Figure  2  was  taken  partly  from  Ref.  2,  partly  from  more  recent 
information  on  OPO's  provided  by  Spectra  Physics,  and  partly  from 
estimates  of  FEL  capability.)  Clearly  in  the  infrared  region  from  3- 
13  pm,  free  electron  lasers  constitute  a  very  large  breakthrough.  In 
between  the  millimeter  waves  and  the  13  pm  infrared^  the 
atmosphere  is  for  all  practical  purposes  opaque.  Naturally  to  be  a 
useful  source  for  atmospheric  sensing,  each  source  would  have  to  be 
fielded,  perhaps  on  a  car  or  truck,  train,  ship,  aircraft  or  space  craft. 
Here  we  briefly  review  some  of  the  important  characteristics  of  each 
source  and  discuss  the  prospects  for  fielding  it. 

1.  Gyro  Devices  at  around  94  GHz 

Gyrotrons  have  been  developed  all  over  the  world^-n.  as  they 
were  developed  for  heating  of  fusion  plasmas  at  frequencies  ranging 
from  35-140  GHz,  they  are  inherently  long  pulse  and  have  high 
average  power  capability.  In  fact,  their  average  power  is  larger  than 
what  is  usually  required  for  a  radar.  Varian  gyrotrons^ 2, 13  are 
available  commercially  with  cw  power  of  300  kW  at  28  GHz,  200  kW 
at  60  GHz  and  100  kW  at  140  GHz.  The  efficiencies  of  the  Varian 
gyrotrons  are  all  around  30%.  A  tabie^^  of  several  Russian  gyrotrons 
with  Gaussian  mode  outputs,  is  shown  in  Table  1.  Also  commercially 
available  gyrotrons  are  available  in  Russia,  from  the  collaboration  of 
the  Institute  of  Applied  Physics  (lAP)  in  Nizhny  Novgorod  and  the 


2 


Tory  Company  in  Moscow.  One  available  there,  which  we  will 
nwntion  shortly  is  the  83  GHz  cw  gyrotron  at  a  power  of  10  kW.  A 
millimeter  wave  radar  with  these  peak  sorts  of  peak  powers,  and  a 
duty  factor  of  0.1-10%,  constitutes  a  breakthrough  in  millimeter 
wave  radar  power  by  several  orders  of  magnitude.  However  it  is 
fairly  modest  as  regards  gyrotron  capability. 

One  of  the  few  gyrotrons  designed  at  a  frequency  optimized  for 
atmospheric  propagation  is  the  NRL  94  GHz  gyrotron  15, 16,  This  has 
generated  up  to  200  kW  of  power  at  94  GHz  with  efficiencies  of 
about  20%.  This  frequency  is  the  clearest  part  of  the  atmospheric 
propagation  window.  Also  tunable  gyrotrons  have  been  developed  at 
NRL.  11 

The  development  of  powerful  gyro  amplifiers  at  frequencies 
near  94  GHz  is  also  procee^ng.  For  many  radar  applications,  an 
amplifier  is  preferable  to  a  free  running  oscillator.  This  allows  phase 
and  amplitude  control  of  the  transmitted  beam.  If  an  amplifier  and 
an  oscillator  are  both  available,  and  all  other  things  are  equal,  the 
amplifier  is  the  preferred  rf  source  to  power  a  radar.  For  one  thing, 
data  processing  is  simpler  if  the  amplitude  and  phase  are  precisely 
known  pulse  to  pulse.  Nevertheless  the  data  processing  can  be  done 
with  an  oscillator  transmitter.  Also,  it  may  be  necessary  to 
coherently  combine  power  from  one  or  more  sources.  This  can  only 
be  done  with  totally  coherent  signals  produced  by  an  amplifier. 
Counterbalancing  the  advantages  of  amplifiers  are  the  fact  that  they 
are  more  complicated  and  expensive  than  oscillators,  due  to  the 
necessity  of  having  the  additional  microwave  plumbing  for  the  input 
signal.  Also  there  are  more  stringent  requirements  of  avoiding  both 
oscillation  and  intercavity  cross  talk  anywhere  in  the  amplifier, 
which  is  usually  longer  ^an  an  oscillator.  For  an  atmospheric  sensor 
application,  the  choice  of  an  oscillator  or  amplifier  will  be  made  in 
terms  of  requirements  as  well  as  price  and  availability  of  whatever 
the  source  would  be. 

At  this  point  the  most  advanced  is  a  gyroklystron  developed  at 
the  Institute  of  Applied  Physics  (lAP)  in  Nizhny  Novgorod  in  Russia^^ 
which  has  generated  power  in  excess  of  50  kW  at  94  GHz.  The  NRL 
quasi-optical  gyrotron,  when  run  in  a  gyroklystron  configuration  has 
also  demonstrated  capability  as'both  an  amplifier  and  a  phase  locked 
oscillatorl8,l9.  The  basic  difference  is  that  the  lAP  gyroklystron  uses 
relatively  low  Q  cavities  so  that  its  bandwidth  is  relatively  high, 
some  fraction  of  a  percent.  The  quasi-optical  gyroklystron  uses  a 
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much  higher  Q  cavity,  so  its  bandwidth  is  much  less.  The  two  are 
complementary  in  that  the  former  could  find  application  where 
bandwidth  is  required;  the  latter  where  great  frequency  and  phase 
stability  is  required. 

All  gyro  devices  at  94  GHz  use  super  conducting  magnets, 
although  future  developments  of  high  harmonic  gyrotrons  might 
render  these  unnecessary 20.  However  if  one  can  live  with  the 
constraints  of  a  super  conducting  magnet,  any  of  these  gyrotrons  are 
reasonably  compact,  certainly  compact  enough  to  be  transported  on  a 
small  truck  or  multi-engine  aircraft.  Thus  they  could  be  fielded  as 
radars.  We  now  discuss  the  sources  in  somewhat  more  detail. 


A.  The  NRL  94  GHz  Gyrotron 

The  94  GHz  gyrotron  is  a  TEis  mode  gyrotron  which  has 
generated  150  kW  of  power  in  a  Ips  pulse  at  an  efficiency  of  15-20%. 
With  appropriate  cooling  on  the  collector  and  window,  it  could 
operate  at  an  average  power  of  10  kW.  This  is  in  fact  only  10%  of 
the  average  power  Varian  has  demonstrated  in  their  140  GHz 
gyrotron.  Another  advantage  to  the  mode  selected  is  that  it  can  be 
fairly  easily  converted  to  a  fundamental  TEn  mode  for  propagation. 

A  mode  converter  has  been  developed  by  General  Atomics,  San 
Diego,  Cal,  and  is  currently  in  place  on  our  gyrotron.  Conversion  to 
the  TEii  mode  has  been  demonstrated.  The  measured  far  field  mode 
pattern  of  the  NRL  94  GHz  gyrotron  in  Fig.(3)  shows  that  the  output 
is  in  fundamental  mode. 

Figure  4  is  a  schematic  of  the  gyrotron  with  the  mode 
converter  in  place,  Fig  5  is  a  photograph  of  the  gyrotron  laboratory. 
Figure  2  shows  the  two  stage  mode  converter  (TEi3-^TEi2  followed 
by  TE21->TEii).  The  equipment  is  compact  and  can  easily  be 
transported  on  a  small  truck.  This  capability  to  radiate  average 
powers  of  order  10  KW  at  94  GHz  into  the  atmosphere  represents  a 
very  new  capability. 

Another  important  figure  of  merit  for  the  gyrotron  oscillator  is 
bandwidth,  which  directly  translates  to  range  resolution  in  a  radar. 
Since  the  Q  of  the  NRL  gyrotron  is  600,  this  corresponds  to  a 
relatively  wide  bandwidth.  The  gyrotron  can  be  tuned  by  varying 
the  voltage  during  a  pulse.  Shown  in  Fig.(6)  is  a  plot  of  the 
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frequency  variation  of  the  NRL  94  GHz  gyrotron,  generated  by  a  time 
dependent  simulation,  where  the  voltage  increased  from  60  to  70  kV 
in  about  20  nsec  (the  time  chosen  for  numerical  convenience).  Here 
the  chirp  is  about  200  MHz,  corresponding  to  a  range  resolution  of 
about  a  meter  for  optimal  data  processing  in  the  receiver  circuit. 
Operation  of  the  NRL  94  GHz  gyrotron  at  different  voltages  have 
basically  confirmed  this  calculation. 

B.  The  Tunable  NRL  Quasi-Opdcal  Gyrotron 

The  quasi-optical  gyrotron  is  a  very  different  sort  of  a  gyrotron 
in  an  optic^  resonator.  It  has  been  developed  at  both  NRL21-24  and 
Lausanne25,  The  beam  propagation  and  the  magnetic  field  are  both 
transverse  to  the  axis  of  the  resonator.  The  mode  output  is  taken  by 
diffraction  around  the  mirror.  The  diffraction  is  characterized  by  a 
fractional  loss  v  per  round  trip  of  the  radiation.  Typically  the  NRL 
experiment  runs  with  v's  between  1  and  5%.  In  theory  and  in  fact, 
the  resonator  operates  only  in  the  lowest  order  fundamental 
transverse  mode. 

Figure  7  is  a  photograph  of  the  quasi-optical  gyrotron.  It  is 
about  as  portable  as  the  conventional  gyrotron.  The  maximum 
power  it  has  achieved  so  far  is  600  kW,  but  at  relatively  low  8% 
efficiency.  More  recent  experiments,  at  the  100  kW  power  level 
have  used  a  number  of  schemes  to  increase  the  efficiency.  The 
maximum  efficiency  achieved  so  far  is  about  20%  without  a 
depressed  collector  and  30%  with  one^^’i^.  While  the  efficiency  of 
the  quasi-optical  gyrotron  is  typically  lower  than  that  of  a  cavity 
gyrotron,  this  is  largely  counterbalanced  by  the  fact  that  because  of 
its  cross  bore  geometry,  a  single  stage  depressed  collector  is 
particularly  simple  to  implement. 

One  of  the  most  attractive  features  of  the  quasi-optical 
gyrotron  is  its  tunability.  It  can  be  tuned  in  three  ways:  •  The 
magnetic  field  can  be  varied  over  a  fairly  wide  range.  The  power  is 
relatively  insensitive  to  field.  Figure  Sa  shows  a  plot  of  mode 
frequency  and  power  as  a  function  of  magnetic  field  for  a  series  of 
lower  power  experiments.  The  device  can  tune  over  about  half  an 
octave  from  about  80  to  130  GHz.  •  The  device  can  be  tuned  with 
the  beam  voltage.  Since  the  frequency  scales  as  B/y,  the  frequency 
can  be  raised  by  lowering  the  voltage.  Figure  8b  shows  the  Voltage 
tuning  curve.  Clearly  the  power  is  now  sensitive  to  voltage  if  only 
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because  the  input  power  is  proportional  to  it.  However  this  voltage 
dependence  does  allow  rapid  tuning  over  a  small  frequency  range  • 
The  cavity  can  be  tuned  by  adjusting  the  mirror  separation,  which 
can  be  done  without  breaking  vacuum.  This  allows  for  very  precise 
tuning. 

It  is  also  possible  that  the  quasi-optical  gyrotron  could  operate 
at  the  second  harmonic.  The  Lausanne  experiment  has  already  seen 
evidence  of  harmonic  operation.  If  it  can  run  with  harmonic 
operation,  it  opens  the  possibility  of  a  single  source  covering  the 
entire  spectral  region  from  80  to  260  GHz. 

C.  The  lAP  94  GHz  Gyroklystron 

The  lAP  amplifier  development  program  is  the  most  mature  at 
this  point^^.  It  has  developed  powerful  tubes  at  35  and  94  GHz 
which  have  been  reported  only  recently.  A  94  GHz  gyroklystron  has 
been  developed  at  lAP  and  is  a  four  cavity  interaction  circuit.  The 
cavities  are  all  in  the  TEqi  mode.  Shown  in  Fig.(9)  is  a  schematic  of 
the  interaction  circuit  superimposed  on  the  background  magnetic 
field.  At  the  peak,  the  field  is  37  kG.  In  Fig.  (10)  is  shown  a  plot  of 
its  power,  gain  and  bandwidth.  In  this  case  the  achieved  power  is  60 
kW  and  the  bandwidth  is  about  400  MHz.  For  the  data  shown  here, 
the  efficiency  was  somewhat  over  30%.  While  the  lAP  program  has 
been  in  existence  for  a  long  time  and  is  quite  mature,  only  recently 
has  information  on  it  been  available. 

D.  The  NRL  Gyroklystron  Amplifier 

This  configuration  uses  a  Fabre  Perot  resonator  with  both  the 
electron  beam  and  magnet  field  direction  perpendicular  to  the  axis  of 
the  resonator.  Just  upstream  from  the  main  resonator  is  a  small, 
lower  Q  prebuncbing  resonator.  At  the  currents  and  perpendicular 
velocities  of  the  beam,  this  upstream  resonator  is  always  below  start 
oscillation  threshold's.  19.  The  experiment  was  done  at  85  GHz 
because  that  was  the  frequency  of  the  available  drive  source.  It  was 
an  extended  interaction  oscillator  (ElO)  which  put  out  a  power  of  1 
kW  for  a  2  ps  pulse.  Power  from  the  EIO  was  fed  into  the  upstream 
resonator,  it  prebunched  the  ■  beam,  and  the  amplified  power  was 
extracted  at  the  main  resonator.  A  schematic  of  the  quasi-optical 
gyroklystron  is  shown  in  Fig.(ll). 
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Since  the  main  resonator  had  a  very  high  Q,  the  bandwidth  was 
very  low,  about  3  MHz,  which  is  about  equal  to  q)/Q.  A  plot  of 
ampliHcation  as  a  function  of  cathode  voltage,  as  compared  with  a 
numerical  simulation  is  shown  in  Fig.(12).  In  this  series  of 
experiments,  the  power  is  about  30  kW  and  the  efHciency  is  about 
10%.  By  increasing  the  current  somewhat,  one  can  exceed  the  start 
current  in  the  main  cavity,  but  remain  below  it  in  the  prebunching 
cavity.  In  this  case  one  can  drive  a  phase  locked  oscillator.  Here  the 
oscillator  is  no  longer  free  running,  but  has  its  phase  controlled  by 
the  much  weaker  input  signal.  This  is  an  alternate  realization  of  a 
phase  controlled  source.  When  run  as  a  phase  locked  oscillator,  the 
device  had  a  power  of  60  kW  and  an  efficiency  of  16%.  In  Fig.(13)  is 
shown  the  phase  locking  bandwidth  as  a  function  of  input  to  output 
power  ratio.  The  curve  is  Adler's  classic  theory26  for  the  case  of 
phase  locking  by  direct  injection  of  power  into  the  output  port  of  the 
device.  This  is 


SfL  =  (fo/2Q)[Pi/Pol»'2 

where  6fL  is  the  locking  bandwidth  and  Pi/Po  is  the  injection/output 
power.  Clearly  the  prebunching  allows  one  to  greatly  exceed  the 
Adler  theory.  Also  because  power  is  injected  into  a  prebunching 
cavity  and  not  the  output,  there  is  no  need  for  a  high  power 
millimeter  wave  circulator. 

The  quasi-opdcal  gyroklystron  is  inherently  a  much  narrower 
band  than  the  conventional  cavity  configuration.  Thus  it  would  find 
its  application  where  great  frequency  stability  is  the  most  important 
characteristic.  An  example  might  be  a  radar  whose  main  role  is  to  do 
Doppler  sorting  of  a  large  number  of  low  cross  section  targets  with 
slightly  different  velocities.  Also  it  does  have  additional  advantages 
due  to  its  simplicity,  intercavity  isolation,  and  adaptability.  None  of 
the  cavities  has  any  inherent  frequency  other  than  c/Lm  where  Lm  is 
the  mirror  separation,  so  it  will  have  a  very  large  inherent  tuning 
range.  The  cavities  are  naturally  isolated  from  one  another,  and  at 
no  point  during  the  experiment  was  there  ever  any  evidence  of  cross 
talk.  Furthermore,  while  the  cavity  Q  was  very  high,  this  too  is 
flexible  and  can  be  lowered  by  moving  the  mirrors  closer  together 
and  by  changing  the  radius  of  the  mirror  edge  so  as  to  vary  the 
output  coupling.  Thus  the  bandwidth  and  the  power  can  be 
increased  with  straightforward  scaling. 
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2  Infrared  Free  Electron  Lasers 


Free  electron  lasers  are  currently  large  facilities  which  take  up 
an  entire  building,  often  a  very  large  one.  As  currently  configured, 
they  cannot  be  fielded.  There  are  two  reasons  the  facilities  are  so 
large.  First  of  all,  they  require  a  multi  megavolt  electron  beam;  and 
second,  because  of  this  beam,  they  need  a  tremendous  amount  of 
radiation  shielding.  Thus  there  are  two  requirements  to  field  an  FEL. 
First  one  must  reduce  the  beam  energy  as  much  as  possible;  and 
second,  one  must  recover  as  much  of  the  beam  energy  as  possible 
with  the  use  of  a  depressed  collector.  There  seem  to  be  two  FEL 
configurations  which  could  potentially  evolve  toward  a  device  which 
can  be  fielded. 

The  first  is  the  electrostatic  accelerator  use  in  the  University  of 
California  at  Santa  Barbara  (UCSB)  FEL^^.  in  this  device,  the  beam 
energy  is  already  recovered  with  about  90%  efficiency  in  a  lasing 
system,  and  higher  beam  recovery  can  be  envisioned  with  better 
beam  optics.  The  beam  energy  in  Ref.(27)  varies  up  to  6  MeV.  With 
the  magnetostatic  wiggler,  the  wavelength  of  the  radiation  varies 
from  about  200|im  to  1  mm.  Thus,  while  this  experiment  recovers 
the  beam  energy,  it  is  still  too  large  kid  it  does  not  operate  at  the 
right  wavelength.  The  key  to  making  it  smaller  is  to  use  a  beam 
energy  of  3  MeV.  In  this  case  the  accelerator  is  about  16  feet  long,  4 
feet  in  diameter  and  weighs  about  4000  pounds,  easily  small  enough 
to  be  flelded  in  a  truck,  train,  ship  or  even  aircraft^®.  Figure  14  is  a 
photograph  of  a  3  MeV  electrostatic  accelerator.  In  order  to  radiate 
in  the  3-13  pm  propagation  window,  the  wiggler  wavelength  has  to 
be  greatly  reduced.  One  way  to  do  this  is  with  the  use  of  an  rf 
wiggler.  One  particular  conBguration,  using  a  quasi-optical  gyrotron 
at  the  second  harmonic  to  form  the  wiggler,  has  been  analyzed.  It 
was  concluded  that  a  portable  system,  which  radiates  about  100 
Watts  of  average  power  in  this  wavelength  regime  could  be 

developed29. 

A  second  possible  configuration  utilizes  a  compact  FEL  such  as 
has  been  developed  at  Los  Alamos^®.  This  employs  a  10-15  MeV 
high  gradient  rf  linac  with  the  beam  produced  by  a  photocathode, 
and  a  small  period  microwiggler.  It  would  generate  tens  of  Watts  of 
average  power  in  the  3-13  pm  wavelength  regime  of  interest.  A 
schematic  of  the  LASL  FEL  is  shown  in  Fig.  15.  The  equipment  is  all 
table  top  size,  and  if  this  were  the  only  consideration,  it  could  easily 
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be  fielded.  Unfortunately  it  is  not  so  simple  because  of  the  issue  of 
radiation  shielding.  The  Los  Alamos  system  is  surrounded  by 
concrete  walls  four  feet  thick  and  is  covered  by  a  steel  roof  six  inches 
thick.  However  this  only  shields  against  the  reflected  bremstrahlung; 
the  main  radiation  shield  is  a  below  ground  beam  dump.  To  field 
such  a  device,  one  would  clearly  have  to  collect  the  beam  energy, 
perhaps  with  a  second  rf  linac  downstream  from  the  PEL  and  phased 
so  as  to  collect  the  beam  energy.  Los  Alamos  has  attempted  this  on  a 
larger  PEL  and  has  managed  to  recover  about  2/3  of  the  beam 
energy  in  a  lasing  system^ I.  One  Stanford  PEL  also  recirculated  the 
beam  energy;  however  this  is  in  a  superconducting  system  which  is 
very  large32.  in  any  case,  if  a  means  can  be  found  to  recover  the 
beam  energy  in  the  Los  Alamos  compact  IR  PEL,  this  could  lead  to  a 
fielded  system  also. 

To  summarize,  gyrotrons  already  exist  in  fieldable  systems,  and 
there  is  every  reason  to  believe  that  IR  PEL’S  can  be  developed  to 
this  point  also. 
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III.  Cloud  Radar  Studies  with  High  Power  Gyrotrons 

High  power  millimeter  wave  sources  open  the  possibility  of  a 
variety  of  active  measurements  of  the  atmosphere.  The  wavelength 
is  one  where  very  few  measurements  have  been  done  up  to  now. 
Millimeter  waves  have  the  difficulty  that  they  are  absorbed  in  the 
atmosphere,  either  weakly  or  strongly  depending  on  wavelength  and 
atmospheric  conditions,  principally  humidity  as  shown  in  Fig.  1.  For 

very  humid  conditions  at  94  GHz,  the  loss  is  about  1  dB/km.  This  is 
fairly  lossy,  but  also  not  so  much  that  high  power  could  not  greatly 
enhance  the  range.  For  instance  by  moving  the  average  power  from 
4  Watts  to  4  kW,  the  additional  propagation  distance  is  increased  by 
15  km  (30  km  round  trip)  in  the  worst  conditions.  Furthermore,  this 
loss  is  only  at  sea  level  where  there  is  considerable  water  vapor  in 
the  atmosphere.  For  propagation  at  angles  to  the  horizontal,  the 
situation  improves  quickly.  As  a  rough  rule  of  thumb,  the  water 
vapor  absorption  at  94  GHz  is  in  the  lowest  2  or  3  kilometers  of  the 
atmosphere. 

Millimeter  waves  have  the  advantage  that  they  propagate  well 
through  particulates,  but  scatter  off  them  strongly  enough  to  give 
easily  detectable  signals,  particularly '  for  high  power  radar 
transmitters.  Convention^  radars  at  microwave  frequencies  really 
do  not  scatter  from  particulates  at  all.  Lasers  radiation,  on  the  other 
hand,  scatters  very  strongly  from  particulates,  but  too  strongly. 

Water  clouds  for  instance  are  virtually  opaque.  A  laser  is  basically  a 
line  of  sight  remote  sensor.  Thus  a  millimeter  wave  radar  could 
resolve  multiple  cloud  layers.  A  low  frequency  radar  would  not  see 
them  at  all,  and  a  laser  would  see  only  partially  into  the  lowest  layer. 
In  this  section,  we  focus  mostly  on  clouds,  because  they  are  very 
important  for  meteorology  and  global  wanning.  However  with  a 
sufficiently  powerful  source,  even  a  clear  atmosphere  can,  and  has 
given  rise  to  detectable  return  signals  at  millimeter  wavelength. 

A.  Millimeter  Wave  Cloud  Radars 

One  of  the  most  promising  applications  for  fixed  frequency 
gyrotrons  as  atmospheric  sensors  is  cloud  radars.  These  have 
potential  application  both  as  ground  and  air  based  systems  for  the 
fundamental  studies  of  cloud  dynamics.  Furthermore,  space  based 
94  GHz  radar  systems  could  give  real  time,  world  wide  maps  of  the 
altitudes,  thickness  and  flow  fields  of  clouds.  This  information  is 
difHcult  to  obtain  directly  in  any  other  way.  Clouds  are  extremely 
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important  in  the  earth's  albido.  It  has  recently  been  pointed  out 
that  the  effect  of  clouds  can  dominate  by  a  large  factor,  the  effect  of 
greenhouse  gases  in  global  warming^^.  jhus  cloud  sensors  must  play 
an  essential  role  in  environmental  research. 

There  are  several  features  of  clouds  which  give  rise  to  radar 
returns.  First  there  is  clear  air  turbulence  in  the  cloud.  This  scatters 
preferentially  low  frequency  radiation,  frequencies  of  10  GHz  and 
lower.  However  in  scattering  from  turbulence  one  is  not  really  sure 
whether  the  turbulence  is  in  the  cloud  or  outside  of  it.  Furthermore, 
it  is  difficult  to  directly  correlate  the  turbulence  with  the  cloud 
physics.  Secondly,  there  are  raindrops,  ice  crystals  and  snow  flakes 
in  the  clouds  which  scatter  the  radar  wave.  However  the  majority  of 
clouds,  particularly  in  their  formation  stages  do  have  these.  Finally, 
there  are  the  aerosols  themselves  which  define  the  cloud.  These 
aerosols  have  very  small  diameters,  typically  about  lOpm.  Since  the 
scattering  cross  section  of  the  individual  aerosol  scales  as  where 
X  is  the  radar  waveleof^th  they  can  only  be  observed  at  millimeter 
wavelength. 

There  have  already  been  a  number  of  studies  of  clouds  using 
94  GHz  radars34-38.  jn  Fig.(16)  is  shown  a  94  GHz  radar  image  of  a 
cloud,  taken  from  L’Hermitte's  Ref.(36).  However  these  radars  all 
use  extended  interaction  oscillators  (EIO's)  and  are  very  limited  in 
average  power,  typically  to  about  4  Watts.  This  greatly  limits  the 
capability  of  the  system  as  regards  to,  minimum  reflectivity  of  the 
cloud,  range  in  humid  air,  and  the  amount  of  data  that  can  be 
accumulated  in  the  cloud. 

In  this  section  we  discuss  the  use  of  the  gyrotron  oscillators  or 
amplifiers  for  remote  sensing  of  clouds,  although  any  particulate 
atmospheric  pollutants  could  be  sensed  by  the  same  principle. 

Clouds  are  made  up  of  small  water  droplets  whose  diameter  ranges 
from  about  10  to  40  pm,  and  whose  droplet  density  varies  from 
about  0.1-100  aerosols  cm'^,  depending  on  the  cloud  type. 

B.  Radar  Scatter  from  Clouds 

Since  the  droplet  diameter  in  the  cloud  is  much  less  than  a 
quarter  of  the  wavelength,  the  scattering  is  in  the  Rayleigh  regime 
and  the  cross  section  is  given  by 
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o=  ji5|k2|D6/A.4. 


(1) 


where  D  is  the  aerosol  diameter,  X  is  the  wavelength,  and  K  =  (m^- 
l)/(m2+2)  and  m  is  the  index  of  refraction.  For  water  at  millimeter 
wavelength,  IK^I  is  given  by  about  0.8.  The  strong  increase  of  cross 
section  with  inverse  wavelength  is  a  potential  advantage  of  a 
millimeter  wave  system.  So  is  the  fact  that  a  tightly  focused  beam 
can  be  produced  with  a  relatively  small  antenna.  Since  the  radar 
beam  intersects  the  entire  cloud,  the  radar  cross  section  is  the  cross 
section  of  each  droplet  times  the  number  of  droplets  in  the  beam. 

The  volume  of  the  cloud  examined  is  the  area  of  the  radar  beam 
times  the  range  resolution  distance,  denoted  h.  For  cloud  studies 
then,  it  is  the  radar  reflectance  t],  the  cross  section  per  unit  volume, 
that  is  the  usual  parameter  used.  Here  11  =  no  and  n  is  the  density  of 
droplets  per  unit  volume.  The  radar  range  equation  becomes  Pr  = 
PthAeTiMicR^,  where  Ac  is  the  effective  area  of  the  receiving  antenna. 

Let  us  consider  the  radar  return  from  a  single  scatterer  in  the 
cloud.  Say  the  incident  radar  pulse  is  given  by 

Et(t)  s  J  df  Gt(f)  exp  2jcift  s  et(t)  .cos[2jcfot  +  (j)t(t)].  (2) 

Here  fo  is  assumed  to  be  the  carrier  frequency  and  e  and  ^  are 
assumed  to  vary  very  slowly  in  time  compared  with  fo.  Each  is 
assumed  to  be  known  and  identical  from  pulse  to  pulse.  The  return 
signal  can  be  constructed  as  in  Section  2.  It  is 

er(t)expi^r(t)  =  s  et(t-2R/c)exp  i(l)t(t-2R/c) 

X  expi[-(47rifo/c)(R  +vt)]  (3) 

where  s  is  proportional  to  the  square  root  of  the  cross  section,  R  is 
the  range  and  v  is  the  dr^  ‘'let  velocity.  Often  the  radar  signal  can  be 
analyzed  by  forming  t'  '  latched  filter  response,  defined  as 

Mn  =  J  dt  et(t-tn)er(t)  exp  i[<j>r(t)  -  (|)t(t-tn)]  = 

Jdt  set(t-tn)  et(ti)  exp  -  <|)t(t-tn)  -47cfoRi/c].  (4) 

That  is  the  matched  filter  response  is  the  convolution  of  the  time 
reversed  transmitted  signal  delayed  by  time  to  with  the  received 
signal.  The  matched  filter  has  two  very  desirable  properties,  one 
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that  it  maximizes  the  signal  to  noise  power,  and  second  that  it  is  a 
simple  way  of  converting  the  bandwidth  of  the  transmitted  pulse 
into  range  resolutions^. 

Consider  the  case  of  a  Gaussian  pulse  with  linear  frequency  chirp.  If 
we  take  e(t)  =  eoexp-[t/T]2  and  <))(t)  =  O.Sat^,  then  we  find  for  the  matched 
filter  response 


Mng  =  (V7c/2)eo2T  exp(4jcifoR/c) 

X  exp{-(tn-2R/c)2[l/2T2  +  a2T2/8]}.  (5) 

The  range  resolution  is  determined  by  the  fall  off  of  the  Gaussian.  If 
this  is  dominated  by  the  frequency  chirp,  the  range  resolution  goes 
roughly  as  caT/V8  «  c/Af  where  Af  is  the  frequency  spread  of  the 
transmitted  pulse.  Note  that  the  matched  filter  response  falls  off 
very  quickly  with  range,  so  targets  can  be  separated  in  range  even  if 
large  and  small  targets  are  much  nearer  to  each  other  than  cT. 

If  an  oscillator  is  used  instead  pf  an  amplifier,  there  is  no 
guarantee  that  the  transnoutted  waveform  will  be  identical  from 
pulse  to  pulse.  However  one  can  still  renormalize  the  transnoutted 
pulse  to  a  desired  transmitted  pulse,  for  each  pulse  in  the  pulse  train, 
and  do  effectively  the  same  thing.  Tnis  was  discussed  in  detail  in 
Ref. (40),  which  introduced  the  notion  of  a  desired  pulse,  which  the 
actual  pulse  is  normalized  to.  While  one  can  use  an  oscillator  and 
renormalize  for  each  pulse,  the  price  to  pay  is  that  much  more  data 
processing  is  necessary  and  also  that  bandwidth  is  not  used  as 
effectively. 


C.  Analysis  of  Radar  Returns  from  Clouds 

Here  we  briefly  examine  the  use  of  a  gyrotron  oscillator  or 
amplifier  as  a  transmitter  for  radar  studies  of  clouds.  We  base  the 
processing  scheme  on  the  matched  filter  response,  as  just  discussed. 
The  cloud  consists  of  a  number  of  individual  scatterers  denoted  by 
the  index  i.  Each  scatterer  moves  with  velocity  v+5vi,  where  v  is  the 
average  velocity  in  the  range  cell  and  5vi  is  a  random  velocity  whose 
average  value  is  zero.  Each  velocity  is  sufficiently  small  that  the 
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scatterers  can  be  regarded  as  stationary  during  the  individual  radar 
pulse.  However  the  velocities  will  be  Doppler  resolved  by  looKing  at 
reflections  from  a  long  sequence  of  pulses.  The  return  signal  can  be 
constructed  by  summing  over  the  returns  from  the  individual 
scatterers.  It  is 

er(t)expi^r(t)  =  Zi  Si  et(t-2Ri/c)exp  i0t(t-2Ri/c) 

X  expi[-(4jcifo/c)(Ri  +vt+8vit)]  (6) 

There  is  not  a  single  transmitted  pulse,  but  a  large  number  of 
pulses  separated  by  interpulse  time  t.  Let  p  be  the  index  denoting 
pulse  number  and  P  be  the  total  number  of  pulses  in  the  pulse  train. 
We  now  calculate  the  matched  filter  response  for  the  p^^  pulse  for 
the  n^  range  interval.  Assuming  that  the  pulse,  in  the  case  of  an 
amplifier,  or  the  desired  pulse  ,  in  the  case  of  an  oscillator,  is 
Gaussian,  the  matched  filter  response  is 

Mnp  =  (Vic/2)eo2T2iSiexp-[(4rcifo/c)(Ri+(v+6vi)px)]  (7) 

where  the  index  n  denotes  the  range  cell  and  p  the  pulse.  We 
assumed  each  droplet  is  at  the  center  of  the  range  cell,  and  the 
summation  is  over  the  droplets  in  the  range  cell.  Associated  with  the 
index  p  is  the  Fourrier  transform  variable  2nz/P  where  z  takes  on 
integer  values  from  zero  to  P*l.  Then  one  can  define 

Qii(z)  =  ZpMnpCxp2jcipz/P  (8) 

To  obtain  the  reflectance,  one  forms  the  summation  ZzlQn(z)l^  ^  Qq. 

In  addition  to  the  summation  over  z,  there  are  four  other 
summations  over  i  and  p  and  over  the  analogous  variables  j  and  q  in 
the  complex  conjugate  summations.  The  summation  over  z  is  a 
simple  summation  of  a  geometric  series  and  turns  out  to  be  equal  to 
P  if  p=q  and  zero  otherwise.  We  then  separate  the  summation  over  i 
and  j  into  a  sum  of  two  summations,  one  having  isj  and  one  having 
i;ftj.  The  former  is  a  simple  summation  over  reflectance's  over  the 
individual  droplets  and  has  N  terms  in  it  where  N  is  the  number  of 
droplets.  The  latter  is  the  sum  of  0.5N(N-1)  terms,  each  one  having  a 
random  phase  assuming  the  droplets  are  uncorrelated.  Thus  the 
ensemble  average  of  this  summation  is  zero.  The  problem  is  that  any 
single  pulse  return  is  not  an  ensemble.  The  rms  value  of  a 
summation  of  random  phases  is  about  N,  so  for  each  value  of  p. 
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the  summation  over  the  i's  and  j's  will  have  about  the  same  order  of 
magnitude.  For  the  case  of  a  sequence  of  pulse  returns,  we  can  often 
exploit  this  summation  over  the  sequence  to  give  the  ensemble 
average.  (In  classical  statistical  mechanics,  the  time  average  and 
ensemble  average  are  assumed  to  be  the  same.)  In  this  case,  the 
summation  over  the  p  index  greatly  reduces  the  value  of  the  i^j 
summation.  In  this  summation,  the  sum  over  p  is  also  a  geometric 
sum  which  can  be  done  simply.  The  result  is  that  this  summation  is 
very  small  unless  5vi-8vj<c/4foPt.  Thus  if  the  velocity  spread  of  the 
droplets  is  sufficiently  large,  the  i^^j  sum  is  small  compared  to  the  i=j 
sum.  If  fo  is  94  GHz,  and  Pt  =  0.1  sec,  the  spread  of  droplet  velocities 
only  has  to  be  large  compared  to  1  cm/s,  a  condition  almost  certainly 
satisfied  for  droplets  in  clouds.  Thus  we  find  that 

Qn=  {0.8P2icWT2d4exp(-2aRn)/16X6Rn4}5;iDi6  (9) 

where  we  have  used  the  relation  between  s  and  o  for  the  droplets, 
the  expression  for  cross  section  in  £q.(l),  and  the  fact  that  the 
antenna  gain  is  x2d2/x2,  where  d  is  the  antenna  diameter. 

To  determine  the  average  velocity  of  the  droplets,  one  forms 
the  summation  lzzlQn(z)|2  s  Fq.  An  analogous  calculation  gives  the 
result 


Fn=  (2PfoViix/c]xQn  (10) 

where  we  have  explicitly  denoted  the  dependence  of  velocity  on 
range  cell.  Other  moments  of  the  distribution  of  velocities  can  be 
found  by  taking  other  moments  of  Qn.  Equations  (9)  and  (10)  show 
how  the  reflectance  and  relative  velocity  in  a  range  cell  are 
calculated  from  the  backscattered  radar  signal. 


D.  High  Power  Millimeter  Wave  Atmospheric  Probes 

We  will  briefly  review  some  of  L’Hcrmittc’s  studies34-37  of 
clouds  using  a  low  power  source,  and  then  discuss  ways  in  which  the 
system  could  be  enhanced  at  higher  power.  Then  we  discuss  briefly 
some  initial  results  on  high  power  millimeter  wave  probing  of  the 
clean  atmosphere. 
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L’Hermitte  studied,  among  other  things,  small  cumulus  clouds 
and  intense  thunderstorms,  both  around  Miami,  Florida.  He  does  not 
use  an  amplifier  as  his  microwave  source,  but  uses  an  oscillator.  To 
do  the  coherent  signal  integration  and  Doppler  processing,  he  uses  a 
coherent  on  receive  approach  and  corrects  the  phase  of  the  received 
signal  according  to  the  phase  of  the  transmitted  signaH^  Except  for 
the  random  phase  from  pulse  to  pulse,  the  EIO  has  extremely  good 
signal  properties.  In  fact,  this  was  the  reason  for  Lhermitte's  choice 
of  an  EIO  over  a  magnetron,  which  could  operate  at  somewhat  higher 
average  power. 

For  a  typical  cumulus  cloud,  like  those  studied  by  L'Hermitte, 
the  droplet  diameter  is  about  10pm  and  the  droplet  density  is  about 
100cm*3,  so  the  reflectivity  is  about  3xl0'l2  cm-i.  With  such  small 
reflectance's  and  low  power  of  the  transmitter,  the  observation  of  the 
cloud  is  difficult.  For  Lhermitte's  case,  with  a  range  of  2  km,  a 
transmitted  power  of  1  kW,  a  range  resolution  distance  h  of  6x10^ 
cm,  and  an  effective  aperture  of  4x10^  cm^,  the  received  power  is 
about  10-13  W.  The  noise  power  of  an  ideal  receiver  at  SOO^K  is 
about  lO'l^W  if  the  bandwidth  is  pulse  length  limited  to  2x10^  s*l. 
However  Lhermitte's  receiver  has  an  actual  noise  of  about  5x10*13. 
His  signal  power  is  below  the  noise  power,  and  he  had  to  integrate 
the  signal  over  a  considerable  number  of  pulses  to  get  good  data.  In 
his  case,  he  runs  the  EIO  at  about  10  kHz  pulse  repetition  rate  and 
integrates  for  3  s,  or  about  3x10^  pulses.  He  finds  that  the 
meteorological  conditions  remain  reasonably  constant  for  these  3s. 

One  of  L’Hermitte’ s  other  observations  was  the  reflection  of  the 
radar  wave  from  a  thunderstorm  in  the  Miami,  Florida  area.  The 
range  at  which  he  was  able  to  operate  was  about  20-30  km. 

However  since  the  air  at  ground  level  was  quite  humid,  he  was  only 
able  to  observe  the  tops  of  the  clouds.  For  instance  one  of  his 
observations  is  at  an  angle  of  34°  and  a  range  of  20  km,  so  that  he 
observes  an  altitude  of  11  km.  He  points  out  that  at  a  30°  elevation, 
the  two  way  absorption  was  12  dB,  whereas  at  a  7°  angle,  it  is  50dB. 

Now  let  us  discuss  the  way  these  cloud  measurements  could  be 
enhanced  at  higher  transmitted  power.  For  the  case  of  the 
thunderstorm,  the  additional  power  could  be  directly  used  to 
overcome  absorption  in  the  humid  air.  For  instance  a  4  kW  average 
power  system  would  overcome  an  additional  30  dB  of  atmospheric 
absorption  and  allow  observation  much  lower  in  the  thunderstorm. 
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Now  consider  the  observation  of  the  cumulus  cloud.  A  small 
cloud  passes  directly  overhead  at  an  altitude  of  about  1.5  km;  as  it 
goes  by,  the  reflectivity  and  vertical  velocities  are  measured  by  the 
radar  at  a  series  of  range  cells  65  m  in  depth.  The  cloud  is  overhead 
for  about  2  minutes.  Since  each  vertical  scan  takes  3  s,  there  are 
about  40  vertical  scans  in  the  cloud  and  about  6-8  range  cells.  A 
higher  power  transmitter  could  enhance  these  measurements  in  a 
number  of  ways.  First,  in  doing  the  same  measurement  as 
L'Hermitte  at  higher  power,  the  radar  beam  could  be  scanned  from 
side  to  side  in  the  cloud  so  that  a  three  dimensional  image  of  the 
cloud  could  be  formed.  Second,  the  range  cell  could  be  reduced  to 
perhaps  10  m  (a  10  MHz  bandwidth)  so  that  the  cloud  could  be 
examined  vertically  in  much  greater  detail.  Third,  the  higher  power 
could  be  used  for  greater  range  to  image  clouds  that  are  not  directly 
overhead.  This  not  only  gives  another  view  of  the  cloud,  but  it  also 
gives  another  component  of  velocity.  For  clouds  with  a  horizontal 
range  of  perhaps  10  km,  the  radar  would  give  a  Doppler  velocity  in  a 
direction  that  is  almost  horizontal. 

The  other  very  important  application  for  94  GHz  cloud  radar 
studies  is  the  proposed  space  based  radar.  Here  the  power  required 
is  greater  than  the  4  Watt  system  used  up  to  now.  Finally  David 
Atlas^2«^3  has  recently  pointed  out  a  number  of  fundamental  issues 
with  regard  to  cloud  dynamics  that  can  be  effectively  studied  with 
94  GHz  radars.  Some  of  these  issues  are  enumerated  in  Table  2. 

We  conclude  by  discussing  a  recent  experiment  on  probing  the 
clear  atmosphere  with  a  high  power  84  GHz  gyrotron"^^.  The 
experiment  was  cw  and  the  transmitter  had  a  power  of  10  kW.  The 
receiver  had  a  sensitivity  of  5xl0'l^  W,  and  the  returned  signal  was 
enhanced  by  using  a  10  second  integration  time.  The  transmitted 
beam  was  scanned  in  angle,  and  at  certain  angles,  there  were  strong 
backscattered  signals.  It  is  most  likely  that  these  signals  are  the 
result  of  particulates,  probably  at  the  top  of  the  atmospheric 
inversion  layer,  at  100-300  meters  altitude.  Clear  air  turbulence  was 
a  possible  scatterer,  but  if  so,  the  turbulence  would  have  had  to  be 
extremely  intense  to  give  rise  to  scatter  at  such  small  wavelength. 
Additional  experiments  are  planned  in  the  future  to  further 
investigate  these  phenomena. 
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IV.  Remote  Sensing  of  Atmospheric  Turbulence 

Here  we  discuss  the  potential  use  of  the  quasi-optical  gyrotron 
as  a  remote  sensor  for  atmospheric  turbulence.  High  power 
gyrotrons  appear  to  provide  a  unique  capability  for  sensing  the  inner 
scale  length  of  the  turbulence  in  an  open  atmosphere  at  fairly  long 
range.  This  is  an  important  property,  because  standard  theories 
relate  it  closely  to  the  strength  of  the  turbulence.  Before  proceeding, 
we  briefly  review  some  properties  of  atmospheric  turbulence  as 
summarized  in  Doviak  and  Zmic^^. 

Power  feeding  the  turbulence  is  provided  at  large  scale  length. 
Nonlinear  convection  of  the  fluid  then  generates  shorter  and  shorter 
scale  lengths  until  the  very  short  scale  motions  are  dissipated  by 
viscosity.  This  intermediate  range  of  cascade  is  called  the  inertial 
range.  The  energy  per  unit  mass  of  the  fluid  per  unit  k  is  d£(k)/dk 
and  has  units  m^/s^.  The  power  per  unit  mass  provided  at  small  k  is 
denoted  e,  and  it  has  units  m^/s^.  For  steady  state  turbulence,  e  is 
also  the  power  dissipated  per  unit  mass  at  large  k,  which  is  also 
equal  to  the  power  convected  through  the  inertial  range  in  k  space. 
The  convection  velocity  through  k  space,  dk/dt,  is  related  to  E  and  e 
through  the  relation  e  =  (dk/dt)x(dE(k)/dk).  The  assumption  now  is 
that  E(k)  depends  only  on  k  and  e.  Assuming  a  polynomial  relation, 
there  is  only  one  which  is  dimensionally  correct, 

I 

dE(k)/dk  =  A  e2/3k-5/3,  (11) 

where  A  is  a  dimensionless  constant  Doviak  and  Zmic  use  the  value 
of  E  to  characterize  the  turbulence.  Table  3  gives  their 
characterization.  Generally,  the  spectrum  seems  to  be  borne  out  by 
experiment. 

The  minimum  value  of  k  is  determined  by  the  power  input  at 
large  scale.  The  m^iximiim  value  of  k  is  determined  by  the  viscous 
dissipation.  The  time  for  the  power  in  this  range  of  k  to  dissipate  by 
viscosity  is  given  by  l/vk2,  where  v  is  the  kinematic  viscosity  in 
m2/s.  The  time  the  spectral  energy  remains  in  a  region  of  size  k  is 
given  by  k/(dk/dt).  Equating  these  two  times,  we  find  that  kmax  is 
proportional  to  el^4/v3/4.  Doviak  and  Zmic  give  the  proportionality 
of  1/5,  so 


kmax  =  0.2el/4/v3/4.  (12) 
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At  sea  level,  v  -  10'^  mVs,  so  for  strong  turbulence,  the  minimum 
distance,  2ic/k  is  about  1  cm. 

Scattering  from  clear  air  turbulence  is  caused  by  fluctuations  in 
the  index  of  refraction  which  in  turn  is  caused  mostly  by  fluctuations 
in  water  vapor  content.  To  determine  the  spectrum  of  index  of 
refraction  fluctuations,  the  spectrum  of  a  scalar  quantity  passively 
advected  by  the  turbulent  field  must  be  calculated^^"^*^.  The  basic 
assumption  of  all  workers,  however,  is  that  the  passively  advected 
scalar  has  the  same  spectrum  as  the  energy,  so  ln(k)|2  =s  Cn^  k*5/3^ 
where  Cn^  is  in  units  of  m*2/3.  Doviak  and  Zmik's  characterization 
the  coefficient  Cn^  in  terms  of  the  strength  of  the  turbulence  is  also 
shown  in  Table  3. 

We  now  consider  the  scattering  from  the  turbulence. 

Assuming  that  5n  is  small,  the  Bom  approximation  for  the  case  that 
the  electric  field  is  evaluated  in  the  far  field  (koD^/r<l,  where  D  is 
the  diameter  of  the  scattering  volume),  gives  the  scattered  electric 
field  at  the  point  r  is  given  by 

Es  *  (ko2Ei/4itr)expikoTjd3T'8n(r’)exp  i(ki-ks)*r’  (13) 

where  Ei  is  the  incident  electric  field  in  the  scattering  region,  ki  is  the 
incident,  and  ks  is  the  scattered  wave  vector.  The  integral  is  over  the 
scattering  volume,  which  is  defined  as  the  intersection  of  the 
patterns  of  the  transmitting  and  receiving  antennae.  Assuming  that 
the  scattering  volume  is  very  large  compared  to  the  scale  size  of  the 
index  fluctuation,  we  find 


IEs|2  =  ko4Ei2VCn2|ki-ksl'“/3/8r2  (14) 

where  V  is  the  scattering  volume.  In  performing  the  integrals  and 
manipulations,  we  have  made  use  of  the  fact  that  the  spectral 
density  per  unit  k  space  volume  is  given  by  In(k)|2/4xk2  for  a 
homogeneous,  isotropic  spectrum.  Also,  we  have  made  use  of  the 
fact  that  the  density  of  states  in  k  space  is  W/(2n)^.  If  the  scattering 
angle  is  0,  then  lki-ksl=2kisin6/2. 

For  backscatter  (e=Jt)  and  k=18(X)m-l  (94  GHz),  the  scattering 
wave  vector  is  almost  certainly  larger  than  kmax>  so  there  should  be 
virtually  no  backscatter.  However  there  will  be  scatter  in  other 
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directions.  The  use  of  powerful  gyrotrons  at  94  GHz  (and  possibly  at 
35  GHz  also)  opens  up  the  possibility  of  examining  the  inner  scale 
length  by  oblique  scattering  measurements.  The  inner  scale  length 
could  be  inferred  from  the  angle  at  which  the  scattering  disappears. 

Let  us  consider  a  specific  example  to  show  how  the  availability 
of  94  GHz  gyrotrons  with  over  100  kW  of  power  can  render  possible 
a  scattering  measurement  of  the  inner  scale  length.  The  scattering 
volume  is  then  roughly  LD^  where  L  is  the  length  of  the  intersection 
of  the  two  antenna  patterns  and  we  have  assumed  the  incident  beam 
is  the  more  tightly  focused  of  the  two.  Then  we  relate  the  field 
amplitude  to  the  power  density  and  add  to  Eq.(14)  additional  factors 
accounting  for  the  atmospheric  attenuation.  Then  Eq.(14)  can  be 
rewritten  as 

Pr  =  ko4AcLPtCn2|ki-ksl-ll/3  exp[-a(ri+rs)]/8rs2  (15) 

where  Pr(t)  is  the  incident  (reflected)  power,  Ae  is  the  aperture  of 
the  receiving  antenna,  a  is  the  atmospheric  attenuation,  and  ri  is  the 
range  from  the  incident  antenna  to  the  scattering  region.  Let  us 
choose  a  configuration  in  which  the  pattered  radiation  is  in  the  far 
Held.  Although  this  is  not  absolutely  necessary,  it  does  simplify  both 
the  analysis  and  the  interpretation  of  any  data. 

Take  a  transmitting  antenna  of  D=3  m  which  can  focus 
radiation  down  to  a  spot  2  m,  2.5  km  away,  and  a  range  to  the 
receiver  of  10  km.  For  relative  humidity  50%  or  less,  the  attenuation 
is  10  dB  or  less.  For  94  GHz  radiation  with  a  3  m  receiving  antenna 
and  an  assumed  inner  scale  length  of  3  cm  (lki-ksl=200),  we  find  that 
Pr/Pi»  2x10-12,  assuming  strong  turbulence,  Cn=3xl0-13,  and  also 
L=D/sin6.  If  the  incident  power  is  lO^  W,  the  received  power  is 
about  2x10-12  W,  4  times  the  noise  power  at  300°  K.  Integration 
over  a  small  number  of  pulses  could  further  enhance  SNR.  Consider 
the  alternative  of  using  an  EIO  having  a  power  of  1  kW  and  an 
average  power  of  4  W.  Now  the  signal  is  a  factor  of  25  below  the 
noise  level,  and  this  is  for  fairly  strong  turbulence.  It  is  unlikely  that 
the  SNR  could  be  significantly  increased  by  integrating  over  a  large 
number  of  pulses;  at  the  shortest  scale  length  (3  cm),  the  correlation 
time  of  the  turbulence  is  undoubtedly  very  short  (10-2s  at  3m/s),  so 
the  scattered  pulses  would  not  be  coherent  with  one  another  for  very 
long. 
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A  configuration  for  bistatic  remote  sensing  of  clear  air 
turbulence  is  shown  in  Fig.  17.  A  transmitter  is  on  the  ground  a 
distance  D  from  the  receiver.  It  points  upward  at  a  shallow  angle  0. 
The  receiver  also  focuses  to  a  point  at  which  it  intersects  the 
transmitted  beam  at  an  altitude  h.  By  varying  the  direction  of  the 
transmitted  beam,  the  angle  6  between  the  transmitted  and  received 
beam  varies,  and  the  scattering  can  be  determined  as  a  function  of 
this  angle. 
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V.  Remote  Sensing  of  Trace  Impurities  in  the  Upper  Atmosphere 

We  now  look  into  whether  the  tunability  of  the  quasi-optical 
gyrotron  can  be  used  to  detect  trace  elements.  At  sea  level,  this 
appears  to  be  an  extremely  difHcult  measurement  due  to  the  fact 
that  the  pressure  broadening  of  the  absorption  lines  is  substantial 
(typically  about  2-3  GHz),  and  the  impurity  concentration  is  very 
low.  That  is  the  absorption  rate  at  the  line  center  scales  as  the 
number  density  of  the  impurity  in  question  Ni  divided  by  the 
collision  rate,  which  is  in  turn  proportional  to  the  ambient  density  Na; 
so  that  the  absorption  rate  is  proportional  to  Nj/Na.  At  sea  level,  Na 
is  very  large,  so  the  absorption  rate  is  correspondingly  small. 
Distinguishing  the  impurity  absorption  from  the  water  vapor 
absorption,  which  dominates  it  by  many  orders  of  magnitude,  does 
not  appear  to  be  a  sinq)le  task.  However  for  impurities  in  the  upper 
atmosphere,  the  pressure  broadening  is  much  less,  so  Ni/Na  is  much 
greater.  Also  there  is  much  less  water  vapor  at  high  altitude  to 
interfere  with  the  measurement. 

The  idea  is  then  to  use  a  high  power  tunable  source  to  do  an 
absorption  (and  possibly  also  phase  shift)  measurement  as  a  function 
of  frequency  near  the  absorption  line.'  (One  might  think  that  at  high 
altitude,  the  molecule  would  spontaneously  re-emit  before  it  is 
collisionally  de-excited.  However  the  time  for  spontaneous  re¬ 
emission  for  rotationally  excited  molecules  varies  from  about  10^  to 
10^  s,  so  this  is  not  the  case^^.)  Measurements  like  this  are  currently 
performed  by  radiometry,  either  ground  or  space  based,  whereby 
the  thermal  spectrum  of  the  atmosphere  is  measured  and  related  to 
the  trace  elements.  An  alternate  scheme  is  to  use  a  millimeter  wave 
source  at  high  power  to  selectively  (as  a  function  of  frequency)  heat 
the  upper  for  diagnostic  purposes. 

A.  Trace  Element  Determination  by  Radiometry 

In  radiometry,  the  thermal  spectrum  of  millimeter  wave 
radiation  is  detected  and  related  to  the  concentration  of  the 
element.^S  These  measurements  can  be  absorption  measurements, 
for  instance  measuring  the  absorption  of  sunlight  at  frequencies  near 
the  absorption  line.  Also  the  thermal  emission  can  be  studied,  and 
this  is  what  we  concentrate  on  here.  For  the  atmosphere  which 
radiates  as  a  black  body  (if  it  were  optically  thick)  the  radiation 
temperature  is  related  to  the  intensity  of  radiation  by  the  Rayleigh- 
Jeans,  or  actually  the  Planck  law.  If  the  radiation  is  emitted  from  a 


22 


region  with  temperature  T  and  s  is  the  path  length,  then  the 
radiation  temperature  Tr  at  the  receiving  antenna  is  given  by 
••  s 

Tr  =  ih  T(s)a(s)exp-T(s),  x(s)  =  iha(s’)ds'  (16) 

where  a(s)  is  the  sum  of  the  absorption  coefficient  of  all  gases  to  be 
considered  and  x  is  the  opacity.  For  a  propagation  medium  which  is 
optically  thin,  as  is  the  case  for  millimeter  waves  in  the  upper 
atmosphere,  the  presence  of  the  gas  of  interest  gives  rise  to  a  slight 
increase  in  radiation  temperature  at  the  frequencies  near  the 
absorption  line. 

The  upper  atmosphere  is  typically  at  a  temperature  of  between 
about  200-3()0oK.  Thus  for  an  optically  thin  propagation  medium, 
the  ratio  of  antenna  temperature  due  to  the  impurity  to  the  source 
temperature  (that  is  the  upper  atmosphere)  is  roughly  the  opacity. 
Ozone  is  a  particularly  attractive  element  to  detect  at  millimeter 
wavelength49,50.  Fig.  18,  from  Waters,  shows  the  zenith  opacity  of 
the  atmosphere  as  a  function  of  frequency,  with  the  absorption  lines 
of  ozone  shown.  Clearly,  there  are  many  possible  lines  to  operate  in 
the  millimeter  wave  band.  For  ozone  at  110.8  GHz,  in  a  standard 
atmosphere,  the  zenith  opacity  is  about  6%,  so  the  antenna 
temperature  due  to  ozone,  at  this  frequency,  is  about  lO^  K  at  the 
line  center.  If  the  line  shape  is  Lorentzian,  as  is  characteristic  of  a 
pressure  broadened  line,  the  temperature  as  a  function  of  frequency 
reflects  this  line  shape.  The  line  width  is  a  function  of  altitude.  If 
sufficient  resolution  of  temperature  as  a  function  of  frequency  is 
available,  the  proflle  of  ozone  density  as  a  function  of  altitude  can 
often  be  inverted.^* 

This  temperature  due  to  the  impurity  is  generally  small 
compared  to  the  atmosphere  temperature.  In  radiometry,  the 
relative  temperature  fluctuation  through  a  channel  of  bandwidth  B  is 
given  by  (BY)' 1^2  where  Y  is  the  integration  time,  so  small 
temperature  changes  due  to  an  impurity  can  be  resolved  by  long 
integration  time.  Parrish  et  al^2  have  recently  studied  trace  element 
concentrations  in  the  upper  atmosphere.  These  are  all  related  to  the 
ozone  cycle,  and  as  such,  the  altitudes  are  about  40  km,  so  the 
pressure  broadening  is  about  10-20  MHz.  They  found  that  at  277 
GHz,  ozone  had  a  temperature  of  8oK  and  the  integration  time  was  an 
hour  to  achieve  SNR  of  600.  They  also  found  that  CIO  had  a 
temperature  about  0.06®  at  279  GHz.  To  achieve  SNR  of  25  took  36 
hours  over  9  days.  The  weakest  line  they  were  able  to  detect  was 
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the  266  GHz  line  of  HO2.  This  line  had  a  temperature  of  about 
0.015<^K  and  took  55  hours  to  resolve  the  line.  Shown  in  Fig  19a  is 
the  result  for  ozone,  and  in  Fig  19b  is  the  result  for  CIO  from  Parrish 
et  al. 


Space  based  millimeter  wave  radiometry  has  demonstrated  the 
capability  of  measuring  simultaneously  the  ozone  and  CIO  profile, 
world  wide  in  real  time.  Waters^  ^  has  presented  beautiful  images  of 
CIO  masses  moving  around  the  world  and  devouring  ozone.  The  real 
time  images  in  Ref  23  change  on  a  time  scale  faster  than  the  time 
required  for  Parrish’s  measurement.  Thus  space  based  passive 
millimeter  wave  radiometry  appears  to  be  the  only  way  to  get  the 
time  resolution  over  the  time  scales  of  interest.  However  these 
measurements  are  themselves  limited  in  some  ways  also.  To  get  a 
world  wide  picture  of  the  ozone  concentration  requires  an  average 
over  altitude,  and  to  get  a  dependence  on  altitude  requires  an 
average  over  a  large  portion  of  the  earth's  surface  taken  over  many 
orbits.  Also  the  radiometers  on  the  space  craft  are  set  at  4 
frequencies,  so  that  the  trace  impurities  to  be  studied  are  set  at  the 
outset  and  cannot  easily  be  varied,  at  least  those  used  in  Ref.  (53). 
Shown  in  Fig.  (20)  are  horizontal  images  of  CIO  taken  over  9  days, 
taken  from  Waters  et  al.  Figure  21  shows  average  vertical  profiles  of 
CIO  from  Ref.  S3,  where  the  data  is  averaged  over  space  and  time. 

We  now  discuss  the  potential  use  of  a  tunable  quasi-optical 
gyrotron  to  do  an  active  measurement  of  upper  atmosphere  trace 
element  impurities.  There  appear  to  be  two  basic  types  of 
measurements  one  can  do,  incoherent  or  coherent.  In  the  former,  the 
phase  information  is  somewhere  lost  on  the  radiation  path,  and  in 
the  former  it  is  not.  Finally  there  is  the  possibility  of  using  the 
quasi-optical  gyrotron  at  high  average  power  for  selective  heating  of 
the  upper  atmosphere  due  to  power  absorption  on  a  particular 
molecular  rotation  line. 

B.  Trace  Element  Detection  by  Incoherent  Scatter 

We  now  consider  the  use  of  a  high  power  tunable  source  for 
measurements  requiring  one  or  two  way  transmission  and  detection 
of  radiation  at  a  set  of  frequencies  near  the  absorption  line.  The 
most  natural  two  way  transmission  path  is  reflection  from  the  moon, 
which  is  visible  from  any  point  on  earth  about  half  the  time.  We 
specifically  envision  the  use  of  a  quasi-optical  gyrotron,  with 
tunability  from  80-260  GHz.  The  radiation  will  be  transmitted  from 
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a  ground  station  on  earth  and  received  in  a  large  number  of 
receivers. 


The  moon  subtends  an  angle  of  about  10*^  radians  as  seen  from 
the  earth.  We  would  like  to  focus  the  radiation  on  as  small  as 
possible  a  spot  on  the  moon,  so  that  the  range  variation  to  the  target, 
and  thereby  the  temporal  spreading  of  the  reflected  pulse  is 
minimized.  At  150  GHz,  a  2  meter  dish  can  focus  the  radiation  to 
about  10%  of  the  moon’s  diameter;  a  20  meter  diameter  millimeter 
wave  qualified  satellite  tracker  can  focus  the  radiation  to  about  1%  of 
the  moon’s  diameter.  Let  us  assume  that  the  radiation  is  reflected 
diffusely  from  the  surface  of  the  moon  into  a  solid  angle  of  about  n 
radians.  Take  for  the  radiated  power,  300kW  and  the  pulse  length, 
10  ps.  Since  the  range  is  about  4x10^  meters,  the  received  power  in 
a  2  meter  dish  is  about  4x10*^^  W  assuming  10  dB  attenuation  on  the 
two  way  path.  If  the  noise  temperature  of  the  receiver  is  about  10^ 
degrees  Kelvin  and  its  bandwidth  is  about  10^  s'^  the  signal  to  noise 
ratio  on  a  single  pulse  is  about  400. 

Thus  the  signal  is  quite  a  bit  above  the  noise.  However  the 
moon  is  a  rough  surface.  That  is  the  surface  scattered  radiation  is 
described  by  for  instance  Eq.(6),  where  the  Si  is  now  a  scattering 
element  on  the  surface.  The  problem  is  that  as  the  frequency 
changes  from  say  the  center  to  the  wings  of  the  absorption  line,  these 
scattering  elements  change  phase  in  an  essentially  random  manner 
due  to  the  fact  that  as  the  frequency  varies  through  the  absorption 
line  (at  least  a  few  tens  of  MHz),  the  phase  of  the  scatterers  are 
randomized.  Thus  in  any  two  measurements,  there  are  order  unity 
variations  between  the  two  received  signals  if  the  measurements  are 
uncorrelated.  The  way  to  resolve  the  absorption  line  then  is  to  do 
an  ensemble  of  measurements.  If  each  measurement  is  uncorrelated 
with  every  other  one,  the  relative  error  for  N  measurements  is 
roughly  jo  decorrelate  the  measurements,  one  would  use 

different  reflection  points  on  the  moon.  (If  the  measurements  are 
partially  correlated  with  one  another,  the  fall  off  in  relative  error 
with  N  is  slower  than  N*  1^2)54,  por  uncorrelated  events,  the 
statistical  properties  are  basically  the  same  as  for  radiometry;  this 
scheme  might  be  considered  radiometry  at  large  signal  strength. 

Thus  for  the  case  of  ozone,  where  the  two  way  zenith  opacity  is 
about  12%  at  110.8  GHz,  roughly  10^  pulses  would  be  required  to 
give  signal  to  noise  ratio  of  10  from  the  center  to  wing  of  the 
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absorption  line.  For  a  single  frequency,  this  would  take  about  a 
second.  To  resolve  the  line  into  ten  frequency  bands  would  take 
about  10  seconds.  For  the  case  of  CIO,  the  two  way  opacity  at  278.63 
GHz  is  about  6x10*^,  assuming  an  ambient  temperature  of  200<^K.  In 
order  to  get  a  signal  to  noise  ratio  of  lO'"^,  about  10^  shots  are 
needed.  If  the  rep  rate  is  10^,  f  s  278.63  GHz,  and  the  parameters 
are  otherwise  as  above,  then  about  3  hours  are  needed  for  each 
frequency  on  the  line.  The  times  here  are  roughly  the  same  as  for 
passive  radiometry,  but  a  much  less  sensitive  receiver  is  needed. 

Now  consider  the  receivers  on  the  ground.  To  achieve  world 
wide  coverage,  there  would  have  to  be  one  transmitter,  most  likely 
on  a  mountain  top  so  as  to  reduce  the  water  absorption  of  the 
transmitted  radiation  on  the  way  up,  and  a  large  number  of  receivers 
on  the  ground.  Wherever  possible  the  receivers  would  be  placed  at 
as  high  an  altitude  as  possible,  and  receivers  could  even  be  air  based 
for  certain  measurements.  Since  the  quasi-optical  gyrotron  is 
tunable  from  about  80-260  GHz,  a  large  number  of  different 
impurities  could  be  measured  simply  by  tuning  the  source.  Thus  the 
recent  availability  of  powerful  millimeter  wave  sources  open  the 
possibility  of  ground  based,  simultaneous,  world  wide  measurements 
of  upper  atmosphere  impurities.  This  is  illustrated  in  Fig.  22a. 


C.  Trace  Element  Detection  by  Coherent  Absorption  Measurements 

We  now  consider  coherent  scattering.  This  can  be 
accomplished  first  by  one  way  propagation,  that  is  the  tunable 
millimeter  wave  source  on  a  satellite  and  the  receivers  on  the 
ground,  or  vice  versa.  Alternatively  it  can  be  accomplished  by  two 
way  propagation.  Here,  millimeter  waves  are  reflected  from  a 
satellite,  and  the  absorption  and  phase  shift  as  a  function  of 
frequency  across  the  absorption  line  provides  information  regarding 
the  trace  elements  in  question. 

We  now  consider  the  use  of  the  quasi-optical  gyrotron  as  a 
tunable  source  for  active  detection  of  ozone.  We  consider  the  110.8 
GHz  absorption  line;  a  frequency  at  which  the  quasi-optical  gyrotron 
can  easily  be  tuned.  To  provide  a  two  way  path,  we  exploit  reflection 
from  an  existing  satellite.  There  are  many  satellites  and  pieces  of 
space  debris  whose  low  earth  orbits  are  known  to  high  precision. 

As  we  will  see  however,  not  all  are  suitable  as  reflectors;  to  be  a 
viable  reflector,  the  size  of  the  satellite  must  be  quite  small. 
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In  active  probing  of  the  upper  atmosphere,  both  amplitude  and 
phase  are  available  to  analyze.  Here  we  consider  only  the  phase, 
since  the  analysis  of  the  amplitude  is  no  different  in  principle  from 
radiometry.  For  a  Lorentzian  line  shape  centered  at  odq,  the  total 
complex  attenuation  and  phase  shift,  E  for  the  two  way  propagation 
is  given  by 


E  =  XV  { ( V  +  i(o)-0)o))/(v2  +  (g)-o)o)2)  } ,  (17) 

where  t/v  is  the  two  way  opacity  at  the  line  center,  and  v  is  the 
pressure  broadening.  At  the  altitude  in  question,  v/2ic  is  about  10 
MHz.  Thus,  for  ozone,  whose  two  way  zenith  opacity  is  12%,  and 
whose  opacity  at  is  about  20%,  the  total  phase  shift  as  one 
crosses  the  line  could  be  as  high  as  100-15°.  If  the  incident  wave  is 
at  frequency  fj  the  phase  of  the  returned  signal  is  given  by 

011”  2jcR(2fi+6f(ii)/c  +  0ti(^“2R/c)  +  A^i(fi) 

+  A4>i(fi+8fdi)  +0  +8d>(fi,p)  (18) 

For  a  satellite  of  known  orbit,  the  Doppler  shift  can  be 
calculated  very  accurately,  and  can  be  as  high  as  a  few  MHz  for  100 
GHz  incident  radiation.  In  Eq.(18),  A<bi(fi)  is  the  ozone  generated 
phase  shift  at  frequency  fi  on  the  way  up,  A<^i(fi+5fdi)  is  the  ozone 
generated  phase  shift  of  the  Doppler  shifted  wave  on  the  way  down, 
0  is  the  phase  shift  from  all  other  components,  determined  from  for 
instance  a  standard  atmosphere  model  and  is  assumed  to  vary  only 
slightly  across  the  absorption  line,  and  5d>(fi,p)  is  the  phase  shift 
generated  by  the  satellite,  where  p  denotes  pulse  number.  The 
satellite  generated  phase  shift  varies  with  frequency  because 
different  frequencies  have  different  scattering  centers;  it  varies  with 
p  because  the  satellite  might  change  its  orientation,  due  to  rotation, 
during  the  pulse  train.  The  difficulty  is  that  5<t>(fi,p)  is  basically 
unknown.  One  can  only  perform  the  measurement  if  5d>(fi,p)  is 
independent  of  both  fi  and  p. 

The  satellite  will  be  considered  to  be  a  set  of  discrete 
scatterers,  with  index  m,  centered  at  R  and  having  position  SRm  from 
the  center.  It  is  the  SRm's  which  contribute  the  essentially  random 
nature  of  the  phase.  If  the  frequency  difference  between  the  waves 


27 


in  the  spectrum  is  denoted  5f,  then  the  contribution  to  the  phase  shift 
is  47c5R5f/c,  where  for  convenience,  we  have  deleted  subscripts.  As 
long  as  this  is  small  compared  to  the  phase  shift  we  are  attempting  to 
measure,  it  can  be  neglected.  This  then  puts  a  limit  on  both  the  size 
of  the  satellite  and  also  on  the  frequency  spread.  For  ozone  at  40  km 
altitude,  assume  a  width  of  20  MHz.  At  a  45°  inclination  angle,  the 
phase  shift  we  are  measuring  is  about  0.2  radians.  This  means  that 
the  satellite  or  debris  radius  must  be  small  compared  to  about  30  cm. 
Also,  the  phase  might  change  pulse  to  pulse  because  of  the  rotation 
of  the  satellite,  so  the  total  measurement  time  must  be  small 
compared  to  a  rotation  period  of  the  satellite,  typically  tens  of 
seconds.  Thus,  a  small  existing  satellite  or  piece  of  space  debris, 
which  is  not  rotating  violently,  could  serve  as  a  reflector  for  rapid, 
ground  based  ozone  measurement. 

We  now  consider  what  it  would  take  to  do  an  ozone 
measurement  at  a  satellite  tracker.  If  we  consider  the  antenna  to  be 
20  m  (a  gain  of  87  dB),  a  bandwidth  which  is  pulse  time  limited  (10^ 
s*0,  a  transmitted  power  of  1  MW,  a  range  of  300  km,  a  5  cm  radius 
spherical  scatterer,  and  a  receiver  temperature  of  lOOO^^K,  we  Bnd 
that  SNR  is  in  excess  of  10^.  This  corresponds  to  a  phase  resolution 
in  a  single  pulse  of  about  (SNR)-i^2.  •  if  20  pulses  within  a  time  of  a 
tenth  of  a  second  to  a  second  are  transmitted  within  the  20  MHz 
interrogation  bandwidth,  the  phase  shift  within  the  absorption  line 
can  be  resolved  and  the  ozone  density  can  be  inferred. 

We  now  consider  a  method  of  making  the  measurement  much 
more  accurate.  The  satellite  generated  phase  shift  can  be  greatly 
reduced  or  eliminated  by  the  use  of  a  specially  designed  satellite. 
Specifically,  we  consider  the  use  of  a  spherical  satellite,  whose  cross 
section  is  independent  of  orientation,  and  whose  cross  section  as  a 
function  of  frequency  can  be  accurately  calculated  and  measured. 

For  trace  elements  such  as  CIO  or  HO2,  discussed  in  Ref.  52,  the  total 
phase  shift  across  the  absorption  line  is  10'^- 10-^  radians.  Thus  the 
satellite  must  be  smooth  to  about  this  fraction  of  a  wavelength. 
However  since  we  are  considering  millimeter  waves,  this  fraction  of  a 
wavelength  means  an  optical  qjality  reflector,  a  standard  capability. 
We  consider  specifically  a  20  cm  radius  spherical  reflector.  In  one 
rocket  launch,  a  number  of  these,  perhaps  ten  or  so,  could  be  placed 
into  orbit  to  provide  nearly  continual,  permanent  coverage.  A 
network  of  satellite  trackers  around  the  earth  would  then  be  able  to 
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provide  nearly  continuous,  world  wide  coverage  of  any  element  with 
an  absorption  line  in  the  millimeter  wave  band. 

It  is  interesting  to  note  that  the  manufacture  of  metal  spheres 
with  this  accuracy  and  launching  them  into  orbit  is  an  established 
technique.  Such  metal  spheres  are  made  by  for  instance  the  Salem 
Specialty  Ball  Corporation  in  West  Simsbury,  Connecticut^^.  They 
have  been  launched  into  low  earth  orbit  as  calibration  targets  for 
measuring  radar  cross  section  of  space  debris  in  the  Orbit  Debris 
Radar  Calibrated  Sphere  project^^  in  1992  and  1993.  (The  micron 
tolerances  in  radius  are  required  so  that  they  can  be  shot  into  space 
from  a  shuttle  based  explosive  launcher.)  The  cost  of  launching  these 
balls  into  a  sufficiently  high  earth  orbit  is  not  that  great,  as  costs  of 
space  flights  go.  Launching  them  would  cost  in  the  neighborhood  of 
one  million  dollars. 

We  now  consider  reflection  from  one  of  these  specially 
designed  satellites  to  measure  trace  impurity  density.  For  instance 
the  242  GHz  CIO  absorptions^  line  would  be  observable  with  the 
quasi-optical  gyrotron  at  the  second  harmonic.  It  has  a  limb 
temperature  of  few  degrees  Kelvin,  and  a  ground  temperature  of 
about  0.05<>K  looking  near  the  horizontal.  Since  an  active 
measurement  has  double  the  path  length,  we  assume  that  at  an  angle 
of  45°  to  the  zenith,  a  value  of  about  O.OS^K,  or  a  phase  shift  of  about 
3x10*4  radians.  The  total  two  way  absorption,  from  sea  level,  at  240 
GHz,  at  an  inclination  of  45°  to  the  zenith  is  less  than  10  dB  for 
relative  humidity  less  than  50%.  If  a  high  altitude  is  chosen  for  the 
measurement,  the  absorption  is  less  still.  We  consider  first  the  use  of 
a  satellite  tracking  station,  with  an  antenna  diameter  of  20  m. 
Consider  a  receiver  temperature  now  of  100°K,  a  total  atmospheric 
attenuation  of  10  dB,  and  a  radiated  power,  assumed  to  be  300  kW  at 
the  second  harmonic.  Then  for  a  10  ps  pulse,  SNR  is  about  6x10^  . 

Ten  pulses  will  then  give  sufficient  SNR  to  accurately  measure  a 
phase  shift  10*4  radians.  To  conclude,  the  tunability  of  the  quasi- 
optical  gyrotron  gives  a  capability  to  rapidly  measure  upper 
atmosphere  ozone  concentrations  with  existing  equipment.  With 
specisdly  designed  satellites  and  accurate  satellite  tracking,  it  could 
give  rapid  measurement  of  trace  impurities  at  much  lower 
concentration.  A  schematic  of  the  measurement  technique  is  shown 
in  Fig.(22b). 

For  what  might  be  considered  the  ultimate  of  such  a  system, 
one  could  place  the  quasi-optical  gyrotron  in  a  high  earth  orbit 
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satellite,  so  that  the  radiation  is  focused  on  the  entire  visible 
hemisphere.  This  would  enormously  increase  the  signal  to  noise 
ratio  since  now  there  is  a  one  way  path.  The  frequency  could  be 
tuned  to  that  of  the  absorber  of  interest  and  the  signal  could  be 
received  at  a  network  of  ground  receiver  stations.  For  instance  if  the 
instantaneous  radiated  power  is  100  kW,  the  received  power  through 
a  2  meter  dish  is  about  10** W,  at  least  6  orders  of  magnitude  above 
the  noise  level  for  a  300° K  receiver  with  megahertz  bandwidth.  If 
the  duty  cycle  of  the  source  is  10*3,  and  the  efficiency  is  30%,  as  has 
already  been  demonstrated,  the  prime  power  needed  for  this 
function  on  the  satellite  is  only  300  W,  assuming  that  it  is  radiating 
all  the  time.  While  the  cost  of  launching  a  high  power  millimeter 
wave  source  into  geosynchronous,  or  other  high  earth  orbit  is  great, 
it  would  provide  a  tremendous  capability  for  rapid,  altitude  resolved 
(by  the  absorption  line  profile),  worldwide  monitoring  virtually  any 
upper  atmospheric  impurity  of  interest.  A  schematic  of  this  scheme 
is  shown  in  Fig.(22c). 

Atmospheric  Modification  With  High  Power  Millimeter  Waves 

We  now  consider  millimeter  wave  sources  at  the  very  highest 
achievable  average  power,  the  several  hundred  kilowatt  to  megawatt 
level,  and  examine  the  possibility  of  atmospheric  modification.  We 
consider  here  the  heating  of  the  atmosphere  by  the  absorption 
(actually  the  differential  absorption)  of  this  radiated  power  by  the 
absorption  line  of  interest.  Let  us  consider  for  instance  the  278  GHz 
line  of  CIO,  for  which  the  opacity  as  measured  by  Parrish  was  about 
3x10*^.  Let  us  imagine  an  average  radiated  power  of  300  kW,  so 
that  the  absorbed  power  by  the  QO  would  be  about  100  Watts.  If 
this  is  radiated  by  a  20  meter,  millimeter  wave  qualified  dish,  and 
the  CIO  which  does  the  absorption  is  spread  out  between  about  30 
and  40  km,  this  then  gives  an  additional  heating  rate  in  the  radiation 
beam  of  about  10*^  degrees  Kelvin  per  second.  If  the  loss  is  due  to 
thermal  conduction  to  the  cooler  gas  outside  of  the  beam,  the 
temperature  rise  of  the  gas  in  the  beam  would  be  limited  to  about 
20K.  As  we  have  seen,  this  temperature  raise  can  be  measured  by 
conventional  radiometry  along  the  boresight.  For  ozone,  a  110.8 
GHz,  the  power  required  to  measurably  heat  the  atmosphere  are 
very  much  less.  Thus  by  using  the  tunable  millimeter  wave  source, 
at  its  highest  average  power  to  heat  the  atmosphere,  one  can  detect 
elements  with  as  low  a  density  as  that  of  CIO. 


30 


VI.  Atmospheric  Sensing  With  Free  Electron  Lasers 

The  field  of  remote  sensing  of  the  atmosphere  with  lasers  is  an 
enormous2*57  ©ne  and  we  will  only  give  a  very  brief  discussion  here. 
As  is  apparent  in  Fig.  (2),  FEL’s,  if  they  could  be  portable  and  fielded, 
could  give  a  breakthrough  capability  in  the  infrared  for  the  case  that 
tunability  is  required.  Tunability  is  required  mostly  for  differential 
absorption  lidar  (DIAL).  Here  the  absorption  of  laser  light  is 
compared  at  the  center  and  far  from  the  absorption  line  of  an 
impurity  in  question.  A  measurement  of  differential  absorption  then 
allows  one  to  infer  the  concentration  of  the  impurity.  There  are 
many  possible  configurations.  First  there  can  1m  a  two  point  path 
measurement  or  a  measurement  at  a  single  point  where  a  remote 
reflector  has  been  placed.  For  either  of  these,  there  has  to  be  a 
cooperative  position  remote  from  the  laser.  Secondly,  the  reflection 
can  be  from  some  typographical  object,  a  building,  tree  or  hillside. 
Since  the  reflection  is  diffuse,  this  second  approach  generally 
requires  more  transmitted  power.  In  either  case,  what  one  ends  up 
with  is  path  integrated  concentration  of  the  impurity.  A  third 
approach  is  to  use  the  existing  atmospheric  aerosols  as  the  scattering 
mechanism.  Here  the  impurity  concentration  can  be  range  resolved, 
but  the  transmitted  power  requirement  is  larger  still.  Furthermore, 
since  the  reflector  is  rough  and  the  propagation  is  affected  by 
atmospheric  turbulence,  the  measurement  is  only  meaningful  when 
averaged  over  a  large  number  of  laser  pulses. 

There  are  many  possible  impurities  which  could  be  detected 
with  IR  DIAL,  but  the  field  is  hampered  by  a  lack  of  tunable, 
powerful  lasers.  To  quote  Fredriksson  58“i>ial  applications  in  the 
mid-IR  region  in  general,  suffer  from  lack  of  tunable,  high  energy 
laser  sources,  and  detection  techniques  is  so  far  not  as  highly 

dedicated  as  in  the  UV  and  visible  DIAL  technique .  A  great  deal 

of  effort  has  been  put  into  IR  DIAL  techniques  by  many  research 
groups,  but  it  is  somewhat  discouraging  that  so  few  Held 
measurements  have  been  made.  However  everything  points  to  a 
breakthrough  when  well  suited  lasers  an  dedicated  detection 
techniques  become  available.” 

In  the  infrared,  the  available  sources  are  basically  line  tunable 
sources.  Thus  DIAL  measurements  have  to  be  the  result  of  a 
coincidence  between  one  of  the  available  laser  lines  and  the 
absorption  line  of  the  impurity  in  question.  For  a  CO2  laser,  the  lines 
are  separated  by  about  50  GHz,  while  the  atmospheric  line  widths 
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are  typically  2-3  GHz.  Thus,  most  of  the  spectral  region  between  9 
and  11  pm  is  uncovered.  For  the  case  where  there  is  a  coincidence  of 
lines,  one  can  do  DIAL  measurements.  For  instance  Killinger  an 
Menyuk59  used  a  CO2  laser  DIAL  system,  where  the  reflected  beam 
was  from  a  typographic  object,  to  detect  CH4.  Shown  in  Fig.  (23), 
taken  from  Ref.(59)  are  the  absorption  as  a  function  of  inverse  wave 
length  for  air  with  and  without  40  ppb  of  CH4.  Note  the  coincidence 
of  the  absorption  peak  of  CH4  with  the  P(14)  line  of  the  CO2  laser. 

With  their  DIAL  setup  and  reflections  from  a  typographic  object, 
concentrations  as  low  as  10  parts  per  billion  were  detected  over  both 
roads  and  airplane  runways.  Their  maximum  range  was  about  2.7 
km.  The  use  of  a  fielded  infrared  FEL  would  open  the  capability  of 
IR  DIAL,  but  without  relying  on  coincidences  between  the  particular 
laser  and  the  absorption  line  of  the  impurity  in  question.  In  addition 
its  continuous  tunability  would  allow  not  only  for  measurements  at 
two  wavelengths,  on  and  off  the  absorption  line;  but  would  also  allow 
for  full  resolution  of  the  spectral  line  shape. 

In  addition  to  its  tunability,  the  FEL  also  brings  an  increase  in 
radiated  power  by  many  orders  of  magnitude,  as  is  apparent  from 
Fig.(2).  This  is  especially  true  between  about  3  and  9  pm 
wavelength,  and  over  a  larger  range  still  if  one  considers 
continuously  tunable  power.  In  DIAL,  the  maximum  range  scales  as 
the  transmitted  power  squared,  since  all  of  the  radiated  beam 
intersects  the  target.  Thus  FEL’s  have  the  capability  of  greatly 
increasing  the  range  of  DIAL  systems.  These  ranges  are  currently 
rather  short,  hundreds  of  meters  to  a  few  kilometers,  and  integration 
times  are  correspondingly  long.  For  instance  Weitkamp^O  et  al 
measured  HCl  concentration  as  a  function  of  distance  from  an 
incineration  ship  by  using  a  low  power  DF  DIAL  laser  system.  His 
results,  taken  from  Ref.(60)  are  shown  in  Fig.  (24).  The  integration 
time  for  each  plot  was  about  an  hour,  and  the  maximum  range  for 
the  system  was  about  2  km.  Thus  if  FEL’s  with  average  power  of 
order  100  Watts  could  be  frelded,  the  entire  area  of  IR  DIAL  could  be 
ripe  for  a  breakthrough. 
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VII.  Conclusion 


Gyrotrons  and  I^L’s  have  potential  for  greatly  increasing  the 
capability  of  atmospheric  sensors.  They  could  extend  weather  radars 
and  lidars  to  new  regimes  of  power,  wavelength  and  tunability. 
Gyrotron  oscillators  at  94  GHz,  and  quasi-optical  gyrotron  oscillators 
or  amplifiers  could  be  fielded  today.  It  is  possible  that  more 
conventional  gyrotron  amplifiers  at  this  frequency  will  soon  be 
available  also.  The  development  of  FEL’s  has  not  yet  matured  to  the 
point  where  they  could  be  fielded,  however  there  are  several 
development  strategies  that  could  lead  in  this  direction. 
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Table  1  —  Capabilities  of  recent  Russian  gyrotrons  at  about  94  GHz. 


RUSSIAN  GYROTRONS  AT  ABOUT  94  GHZ 
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FLOW  IN  YOUNG  CONVECTIVE  CLOUDS;  EFFECT  OF  TURBULENCE 
ON  CONDENSATION 

MIXING  ACROSS  INVERSION  LAYERS 


Table  3  —  Characterization  of  the  strength  of  clear  air  turbulence. 


Turbulence 

Cn2  m-2/3 

e  m2/s-3 

Weak 

6x10*17 

3x10*3 

Moderate 

2x10*15 

8.5x10*3 

Severe 

3x10*13 

6.8x10*2 
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Fig.  1 
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~  Millimeter  wave  atmospheric  attenuation  and  phase  shift  as  a 
function  of  a),  frequency  and  b),  humidity,  from  Ref.  1 


M  ■ 


(%)  NOtSSIVVSNVUl  (r)  AOU3N3  iSlOd  U3SV1 


Fig.  2  —  a)  Atmospheric  attenuation  in  the  infrared  as  a  function  of 
frequency,  b)  Approximate  capability  of  different  laser  systems. 
Lines  ending  in  dots  are  line  tunable,  lines  ending  in  stars  are 
continuously  tunable. 
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Fig.  3  —  Far  field  mode  pattern  of  the 


94  GHz  GYROTRON  SCHEMATIC 


45 


Fig.  4  -  Schematic  of  the  NRL  94  GHz  TEn  gyrotron 
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Gun  Voltage  (kV) 

Fig*  S  Tuning  of  the  NRL  quasi-optical  gyrotron  with  a)  magnetic  field, 
and  b)  voltage 
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Fig. 


10—  Power,  gain  and  frequency  response  for  the  lAP  94  GHz 
gyroklystron 
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DEPRESSED 


quasi-optical  gyroklystron 
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□  Simulation 
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Cathode  Voltage  (kV) 

Fig.  12  —  Gain  as  a  function  of  cathode  voltage  for  the  NRL  quasi-opdcal 
gyroklystron  amplifier 
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gyroklystron  phase  locked  oscillator 


55 


;itu’ 


57 


Fig.  17  -  Configuration  for  ground  based  remote  sensing  of  clear  air 

turbulence  from  ground  level  up  to  an  altitude  of  several  kilometers. 
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18—  Zenith  opacity  of  the  atmosphere  as  a  function  of  frequency  with 
and  without  the  water  absorption*  (without  water  corresponds  to 
very  high  altitude).  The  lines  are  the  absorption  due  to  ozone. 


0.0  ar~i!o^i^^^o  C'O  (ppb.v.) 

20  —  Variation  in  .e  distribution  of  QO,  from  Waters  et  al,  during 
January  2-13,  1992  CIO  abundance  is  indicated  by  the  grey  scale 
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Rg.  22—  Schematic  of  a  variety  of  active  upper  atmosphere  trace  element 
detection  schemes,  a)  reflection  from  the  moon,  b)  reflection  from  a 
spherical  callibration  satellite,  and  c)  a  tunable  millimeter  wave 
source  in  space. 
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